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The criterion for the maximum incoherent conductivity a c ,max in the anisotropic system is shown to be 
o'cmaxO'ab = C (e 2 / h) 2 U2d , where U2D is the two-dimensional electron density, and C is a numerical constant 
of order one. A model of c-axis interlayer conductivity of layered system in the incoherent regime is developed. 
The temperature dependence of the c-axis conductivity is derived to be o~ c (T) ~ T a , where a is determined 
by the effective strength of the interplane voltage fluctuations. The "pseudogap" phenomenon of suppressed 
interplane conductivity at low frequencies (or low voltages) is discussed. 

PACS numbers:74.25.-q, 74.25.Fy, 74.25.Gz 



The nature of the interplane transput properties 
of various layered materials (cupratesoo, ruthenate 
Sr2RuO/(ti, n and some other inorganic and organic 
compounds^) has been of great experimental and the- 
oretical interest. In these systems, the incoherent c-axis 
transport (defined as non-metallic temperature depen- 
dence of resistivity and the non-Drude optical conduc- 
tivity) can coexist with metallic in-plane transport. For 
most of the cuprate compounds the c-axis conductiv- 
ity is incoherent, although for some cuprates (especially 
upon ovprdoping) the c-axis conductivity is metallic or 
coherento. Another example is Sr-2Ru04, where the in- 
coherent c-axis conductivity is observed above 130 Ku. 

The interplane tunneling must be viewed as a many- 
body process from a class of orthogonality catastrophes 
when the tunneling electron unavoidably excites various 
degrees of freedom of the conflated electron liquid. In 
this paper I develop a modeH to describe the incoher- 
ent interplane tunneling which takes into account the 
effect of charge fluctuations excited during the process 
of tunneling. The effect of these charge fluctuations is 
to decohere and to suppress strongly the probability of 
interplane tunneling. 

The coherent band c-axis transport (for instance, the 
case for Sr2RuOi below 130 K) is well described and un- 
derstood in the framework of anisotropic Drude model, 
while the theory of, the incoherent c-axis transport re- 
mains controversiali'BI. Without relying on any specific 
description of the incoherent tunneling, it is possible to 
derive a simple estimate for a crossover between coherent 
and incoherent regimes. This estimate can be made by 
comparing the duration of the interplane hopping time 
Tf- lop and the in-plane scattering time T a b- In the case 
of incoherent electron hopping the c-axis diffusion con- 
stant D c can be written as follows D c = dr/rhop, where 
d is the interplane distance. The Einstein relation allows 
to relate the c-axis conductivity and the c-axis diffusion 
constant 

O-c = ~j = , (1) 

where V2D — to*/ (^h 2 ) is the two-dimensional density of 



states. The in-plane scattering time T a b can be estimated 
from the standard Drude formula a a b — {e 2 V2DD a b) / d, 
where D a b = VpT a b/2- Thus the hopping time Thop 
and the in-plane scattering time T a b can be calculated 
from experimental values of the corresponding conduc- 
tivities (and the density of states v^d)- The descrip- 
tion of incoherent interplane hopping is consistent if the 
hopping time Thop is much longer than the scattering 
time Tab. It exactly means that the interplane tun- 
neling is "disrupted" by the inelastic scattering in the 
planes. Due to the long duration of the interplane tun- 
neling (Thop 3> Tab), the in-plane scattering influences not 
only the in-plane propagation but also independently the 
tunneling process. In other words, collective excitations 
(e.g. voltage fluctuations) in a wide range of frequen- 
cies larger than h/Thop can be excited by the tunneling 
electron. Note that the opposite condition is required for 
the validity of the coherent hopping, namely, the in-plane 
scattering time r a b must be longer than the characteristic 
interplane tunneling time h/t±, where t±_ is the tunneling 
matrix element. Apparently, the crossover between two 
regimes happens when the interplane hopping time Thop 
becomes comparable or equal to the in-plane scattering 
time T a b- By multiplying the left and right hand sides of 
the expressions for the c-axis and in-plane conductivities 
given above and using the crossover condition Thop — i~ a b, 
we derive an estimate for the maximum incoherent c-axis 
conductivity <y c ,max 

(2 \ 2 
-r\ n 2D, (2) 

where n2D = k F /(2ir) is the two-dimensional electron 
density, C is a constant of order 1 which cannot be de- 
termined with the above qualitative argument. The con- 
dition (||) was derived from the premises of the incoherent 
transport. Accordingly, it can be compared with a very 
similar criterion for the minimum metallic conductivity 
derived from the anisotropic Drude model with simpli- 
fying assumptions!! The condition (^) can be consid- 
ered as the generalisation of the Ioffe-Regel-Mott min- 
imum metallic conductivity for anisotropic layered sys- 
tems. The Eqn. can be rewritten as following 
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where cfmir — \kp is the Mott-Ioffe-Regel conductiv- 
ity. The salient point is that the crossover from the in- 
coherent conductivity to the metallic c-axis conductiv- 
ity can happen for a c significantly smaller than (Jmir 
(if a a b 3> tmir), contrary to the naive application of the 
standard Mott-Ioffe-Regel condition. The above fact that 
at crossover a c <C <Jmir was noticed widely for high-T c 
cuprates. 

It is interesting that the condition (|J) is expressed only 
through the two-dimensional (2D) density n 2 D and inde- 
pendent from interaction related renormalisations. For 
a typical metallic 2D density n 2 u ~ 10 15 cm~ 2 of lay- 
ered materials and choosing C ~ 1, the condition for the 
product of in-plane and out-of-plane conductivities gives 

Oc,rnaxOab ~ 10 6 (f2cm)~ 2 (or equfvalently, p c ,minPab ~ 

10 _6 (f2cra) 2 ). This criterion is roughly consistent with 
experimental observations of crossover between incoher- 
ent hopping and coherent band transport. Indeed, the 
crossover as a function of temperature in Sr 2 RuO^ and 
the crossover as a function of doping in various cuprates 
occurs roughly when p c p a b ~ 10 -6 (flcm) 2 . 

Next I demonstrate briefly the generalisation of the 
tunneling formalism in the presence of the V|Olta|ge fluc- 
tuations (more details can be found in RefBQEl). The 
voltage fluctuations which suppress strongly the tunnel- 
ing can be related to the in-plane conductivity <r a b(w) 
by the use of the fluctuation-dissipation theorem. The 
main part of the paper is devoted to the calculation of 
the specific results for the temperature and frequency 
dependencies of the c-axis conductivity <t c (uj, T). At the 
end of the paper, I discuss experimental data which can 
be understood with the help of the developed here theory. 

The time-dependent tunneling formalism allows us to 
derive the expression for the tunneling current I(V) be- 
tween two metallic planes as a function of applied voltage 
V: 

I(V) = — dEdE'dkdk'\t ± (k, k')\ 2 Ai(k, E)A 2 (k' , E') 

x {f(E)(l - f(E'))P(E + eV-E',k- k') - 

- /(£')(! - f(E))P{E' -eV-E,k- k')}, (4) 

where Ai^ 2 (k,E) are the spectral functions, f(E) is a 
Fermi function, and P(E, k) is a Fourier transform of 
the tunneling probability function P(r — r',t — t'). The 
essential new quantity here in comparison with standard 
treatmento is the tunneling probability function P(r — 
r',t-t'): 

P(Sr = r -r',8t = t-t') = exp{-R(Sr, St)),R(Sr, St) = 



dujd 2 q 



< SV^ > (1 - cos(u6t + q6r))coth—, (5) 



where < SV 2 ^ > is the spectrum of the interplane volt- 
age fluctuations. The tunneling probability function 
P(Se, 8k) is the probability of loosing or gaining energy 



Se and momentum 8k from "the electromagnetic envi- 
ronment" of the voltage fluctuations during tunneling. 
Perhaps, the most serious assumption in the derivation 
of Eqn. |] is that the effects of interactions (e.g. voltage 
fluctuations) on the tunneling probability P{8E, 8k) and 
the in-plane spectral functions A\ y2 {k,E) can be consid- 
ered separately. 

Given the spectrum of the interplane voltage fluctu- 
ations < 8V 2 U >, the tunneling probability P(Sr,St) 
(or a Fourier image P(Se,8k)) can be calculated from 
Eqn. ^|. Under certain requirements the tunneling prob- 
ability function P(Se, Sk) is independent of the change of 
the momentum Sk implying the diffusivity of tunneling 
(the momentum of the electron parallel to the planes is 
not conserved after the tunneling). These requirements 
have transparent physical meaning to ensure the diffu- 
sivity of tunneling (see Ref.utf). The tunneling must 
be "slow" Thop 3> h/kT, and the in-plane conductivity 
a 2D = Vabd cannot be much larger than the quantum 
conductance oq = e 2 /h (therefore the spreading of the 
electron density from the tunneling point is also slow). 
The calculations are significantly simplified if the tun- 
neling is diffusive. We consider below only the case of 
the diffusivertHnneling, since these conditions are satis- 
fied (see RefBS) for the most interesting case of cuprates 
. If any of these conditions (still in the regime of incoher- 
ent hopping T^p > T a ii) is violated, the full structure of 
the probability function P(Se, Sk) is necessary to consider 
to account for a partial momentum conservation. 

Under the conditions of the diffusive tunneling, Eqn. ^ 
can be simplified by the integration over the in-plane mo- 
menta (fc, k'): 
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In the derivation of the Equation ^ we assumed addi- 
tionally that: (a) the density of states V2d(E) — v 2 d is 
constant close to the Fermi surface (b) the tunneling ma- 
trix element t± is given by the averaged value around the 
Fermi surface (under conditions of the diffusive tunnel- 
ing, thej-angular variation of the local tunneling matrix 
elementL3 is not particularly important). 

To model the interplane voltage fluctuations spectrum 
< SVgu > we consider a pair of planes coupled by 
Coulomb interaction. It can be shown that the necessary 
generalisation to a stack of planes coupled by Coulomb 
interaction does not change the essential results. The 
voltage fluctuations which most effectively "detune" (or 
decohere) the interplane tunneling are antisymmetric 
with respect to the planes. When electron tunnels be- 
tween planes, these antisymmetric fluctuations are read- 
ily excited and dampen the tunneling. The voltage fluc- 
tuations (Johnson-Nyquist noise) can be calculated by 
using the fluctuation-dissipation theorem which relates 
them with the frequency dependent in-plane conductiv- 
ity a 2D {uj) = a ab (Lu)M: < SV 2 W >= ^U a 2 (g)/m X ^, 
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where x£ and V a {q) = 2ire 2 (\ — e~~ qd )/q are the charge 
susceptibility and Coulomb interaction in the antisym- 
metric channel in two dimensions. The integration over 
momentum gives the local voltage fluctuations 



(jq (2-KReai D (w)q 2 c d\ 

where oq — e 2 [Ti is the quantum conductance^, and q c 
is an upper momentum cutoff which we can take as the 
inverse screening length k — 2 r Ke 2 V2D or the inverse of 
the interplane distance l/d. 

Under the assumptions of the diffusive tunneling it 
is sufficient to calculate the local tunneling probabil- 
ity P(Sr = 0,5t), because the tunneling probability is 
strongly peaked at Sr = 0. The function R(5r = 0, St) 
in the exponent of the P(5r — 0, St) is R(5r — 0, St) ~ 
p- + aln(2n<j2DK 2 d/(kT)), where 
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and = ho~id / '(2<jQkT). The first term in the expres- 
sion for R(Sr = 0, St) is determined by the low frequency 
fluctuations (ui < kT), while the second term is due 
to the higher frequency fluctuations (w > kT). The 
high-frequency cutoff uj c can be taken as either the in- 
verse of the underbarrier traversal time h/T tr ~ leV or 
the frequency of the interband transitions (~ 1 — 2eV) 
above which the c-axis conductivity increases substan- 
tially. The tunneling probability P(Sr — 0, St) allows 
to calculate the temperature dependence of the c-axis 
conductivity from Eqn. |g[ If the low frequency in-plane 
conductivity is not strongly temperature independent 
Cafc(w — *■ 0) = const (e.g. the residual conductivity), 
the essential temperature dependence of the c-axis con- 
ductivity is 



a c {T) ~ T° 



(10) 



where a (defined by Eqn. |9|) is very weakly temperature 
dependent exponent of the integrated weight of the volt- 
age fluctuations. The above result for the temperature 
dependence of the c-axis conductivity as well as the fre- 
quency dependence discussed below are the main results 
of the paper. It is important to notice that the main 
contributions to the exponent a (as can be clearly seen 
from the Eqn. [9]) come not only from small frequencies 
but also from high frequencies where the optical in-plane 
conductivity <Jid{w) is small. 

The results for the temperature dependence of the c- 
axis conductivity depend sensitively on the low frequency 
fluctuations spectrum (or in other words, on the temper- 
ature dependence of the in-plane conductivity a a b{uj < 
kT)). For instance, if as observed for the optimally doped 



cuprates the in-plane conductivity has a strong temper- 
ature dependence a a b ~ 1/kT with a negligible residual 
conductivity, the temperature dependence of the inter- 
plane conductivity is modified to er c (T) ~ T Q_1 . An- 
other complication is that the low frequency spectrum of 
interplane voltage fluctuations are affected by the small 
but finite interplane conduction. Due to the interplane 
conduction, the interplane voltage fluctuations are sup- 
pressed for the frequencies below the inverse hopping 
time (u < 1/r/jop). If the condition l/r^ op <C kT is sat- 
isfied which is the case only for the strongly anisotropic 
materials, we can neglect the effects of the interplane con- 
duction on the fluctuations spectrum. But if the hopping 
frequency 1/rhop approaches kT, then the suppressed in- 
terplane fluctuations below ui < 1/t^ op do not dampen 
the tunneling. The account of self-consistent corrections 
from such effects (together with a possibility of tunneling 
with partial momentum conservation) presents a com- 
plicated problem, and this paper is intended to capture 
only the main results under the simple conditions. For 
the same reason, the description of the crossover from the 
incoherent hopping to the coherent band propagation ap- 
pears to be a difficult problem. 

It is interesting to note that the formalism using the 
tunneling probability function can be used to rederive the 
anomalous corrections-to the tunneling density of states 
of the disordered filmEHl. The voltage fluctuations in the 
disordered film are given by the expression of Eqn. with 
V a (q) substituted by Coulomb interaction in the two- 
dimensional plane V2d(q) = 27r/ g. For the disordered 
film with a short elastic scattering time r the in-plane 
conductivity 02,0(0;) is temperature and frequency inde- 
pendent (up to the frequency 1/r). The calculation of 
the tunneling probability gives 
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This nonperturbative expression for the diffuspp 
anomaly is equivalent to the results in the literatureEj. 
The fluctuation-dissipation theorem used above guaran- 
tees the equivalence of the two approaches, the one based 
on the self-consistent linear response of the electron liq- 
uid to the tunneling electron and the other based on the 
account of the voltage fluctuations. 

One of the most distinct properties of incoherent trans- 
port is the frequency dependence of the optical conduc- 
tivity. The c-axis optical conductivity a c (uj) in the inco- 
herent regime is generally flat (or almost frequency in- 
dependent) with suppression at low frequencies in some 
cases. Even qualitatively, it is very different from metal- 
lic (Drude peak centred around w — 0) or localisation re- 
sponses (a peak at finite frequency with <t(uj — > 0) ~ 0). 
The optical conductivity can be calculated in the same 
framework of the inelastic tunneling formalism. Under 
the conditions of diffusive tunneling, the functional de- 
pendence of the conductance u c (V) = I(V)/V on the 
applied voltage is equivalent to the frequency depen- 
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dence of the optical conductivity a c (ui)B. The c-axis op- 
tical conductivity <t c (lo) can be calculated numerically 
(see Fig. 1). The magnitude of the conductivity <t c (uj) 
is arbitrary in Fig. 1, because the tunneling matrix ele- 
ment varies for different materials (while there is an up- 
per bound on the incoherent conductivity (Eqn. [^), there 
is no lower bound). The suppression of the c-axis con- 
ductivity at lowpfrequencies is due to the soft Coulomb 
blockade effectQEJ, namely, the tunneling is suppressed 
by the slow diffusive spreading of the electron from the 
tunneling point. The frequency scale of the suppression 
of the conductivity is essentially given by 2ir<j2DK- 
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FIG. 1. The c-axis optical conductivity a c {ui) as a function 
of frequency u for two different temperatures. The param- 
eter a is chosen to be 1. The lower curve is for the lower 
temperature. 

Some systematics of the interlayer conducting proper- 
ties of the cuprates can be understood from the view- 
point of the Coulomb blocking of the tunneling by the 
voltage fluctuations. The most appropriate materials 
(not complicated by the proximity to the regime of the 
anisotropic band transport) to test the theory are the 
strongly anisotropic cuprates like Bi — 2201 and Bi — 2212 
(only inter-cell hopping which determines the transport 
properties is considered here). These materials have the 
incoherent c-axis conductivity (clearly seen by the non- 
metallic temperature dependence of the conductivity) 
throughout most of the phase diagram as a function of 
doping. The effective integrated weight of the voltage 
fluctuations a (Eqn. ^) can— be determined from the in- 
plane optical measurements!^. Then the temperature de- 
pendence of the conductivity cr c (T) ~ T a can be tested 
experimentally. The estimates for the parameter a for 
various cuprates give values I < a < 5 (up to the uncer- 
tainty of the upper cutoff frequency u> c ), which are consis- 
tent with the exponents q£ the temperature dependence 
of the c-axis conductivity.!!! Unfortunately, as mentioned 
above it is difficult to apply the theory to some cuprates 
in the proximity of the crossover to the anisotropic band 
propagation (as a function of the doping or temperature), 
because the voltage fluctuation spectrum depends sensi- 
tively on various factors. 



The pseudogap suppression of the c-axis conductiv- 
ity is explained by the proposed model to be due to 
the Coulomb blocking of the tunneling. Because with 
underdoping the in-plane conductivity a a b is reduced, 
the "pseudogap suppression" of a c (u>) is consistently ex- 
pected to be more pronounced due to the increase of 
the parameter a (in agreement with Eqn. ^|). Since 
the c-axis conductivity depends exponentially on the 
parameter a, the c-axis conductivity a c is suppressed 
stronger than the in-plane conductivity upon underdop- 
ing in agreement with experimentally observed increase 
of anisotropy p c / p a b — &ab/&c- The interlayer tunnel- 
ing spectroscopyr-ejxperiments provide a support for the 
proposed pictured. Other aspects of the model are con- 
sistent with optical measurements ofj-the conductivity 
ct c (lu) in the normal state of cupratestJ: in addition to 
the "pseudogap" suppression at low frequencies the c- 
axis conductivity ct c {uj) is suppressed in a wide range 
of frequencies upon lowering the temperature, and the 
spectral weight transferred to frequencies above the bare 
tunneling matrix element frequency t±/%. 

The understanding of a whole set of phenomena asso- 
ciated with the pseudogap observed for the underdoped 
cuprates is essential to describe consistently the c-axis 
transport properties of cuprates, and this problem is not 
addressed here. The applicability of the proposed theory 
to cuprates is quite limited, because the in-plane den- 
sity of states V2D was assumed constant. Nevertheless, 
the c-axis "pseudogap" suppression of the conductivity is 
shown to arise due to the soft Coulomb blockade even for 
the constant in-plane density of states close to the Fermi 
surface. Possibly, it describes the c-axis optical conduc- 
tivity of Sr2RuO^ above 130 K (no apparent pseudogap 
suppression was observed for the in-plane properties in 
this material). 

In conclusion, the incoherent interlayer conductivity of 
layered materials is proposed to be due to the dampening 
of the interlayer tunneling by the voltage fluctuations. 
The power-law for the temperature dependence of the c- 
axis conductivity is expected, with the exponent of the 
power-law determined by the integrated weight of the 
voltage fluctuations. The "pseudogap" suppression of the 
c-axis conductivity at low frequencies (or low voltages) 
results from the Coulomb blocking of the tunneling. 

This work has benefited tremendously from fruitful dis- 
cussions with A.J. Leggett. 
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